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Abstract 

The homogeneous Lippmann-Schwinger integral equation is solved in momentum space by using 
confining potentials. Since the confining potentials are unbounded at large distances, they lead to 
a singularity at small momentum. In order to remove the singularity of the kernel of the integral 
equation, a regularized form of the potentials is used. As an application of the method, the mass 
spectra of heavy quarkonia, mesons consisting from heavy quark and antiquark (Y(bb),ip(cc)), are 
calculated for linear and quadratic confining potentials. The results are in good agreement with 
configuration space and experimental results. 
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I. NUMERICAL SOLUTION OF NON-RELATIVISTIC EQUATION 



The solution of non-relativistic and semi-relativistic Schrodinger equation with confining 
potentials is interesting in various phenomena in physics, from particle to atomic physics. 
Many numerical methods have been developed to study such systems in configuration and 
momentum spaces, such as recent asymptotic exact solution of two- and three-body problems 
l|, l2j . In this work we use a regularization method to study the two-body systems, which 
interact by confining potentials. The bound state of two equal mass particles in momentum 
space and in a PW representation is described as: 

1 f°° 

Hp) = —-i / dp'^V^p')^'), (1) 
V l (p,p') = 2<K J dxP^V&p^x). (2) 

Since the integral equation ([I]) is singular for confining potentials, consequently the calcu- 
lated energy eigenvalues would not be in agreement with the exact analytic binding energies. 
To overcome this problem one can use the regularized form of confining potentials to remove 
the singularity of the kernel. To this aim one can keep the divergent part of the potential 
fixed after exceeding a certain distance, which creates an artificial barrier. The influence 
of tunneling barrier is manifested by significant changes in the energy eigenvalues at small 
distances. For numerical solution of integral equations ([1]) and (2) and for discretization of 
continuous momentum and angle variables we have used Gauss-Legendre quadrature grids 
with hyperbolic plus linear (200 mesh points) and linear (100 mesh points) mapping cor- 
respondingly. The momentum integration interval [0, 00) is covered by a combination of 
hyperbolic and linear mappings of Gauss-Legendre points from the interval [—1, +1] to the 
intervals J0,pi] + \puP^ + ^P^ as: 

hyperbolic linear 

1 + X _ P3~P2 P3 + P2 

-^hyperbolic i , 2 1 \ ? -^linear ^ 2 

pi \ Vl Pi' x 

The used values for p ly p 2 and p 3 in our calculations are 1.0, 3.0, 10.0. In the following 
we present the calculated energy eigenvalues for linear and quadratic harmonic oscillator 
potentials and we investigate the agreement to energy eigenvalues obtained from analyti- 
cal solution of Schrodinger equation. Fourier transformation of regularized form of these 



potentials to momentum space is given by: 



V(r) = air : V(p,p',x) = air c 5 3 (q) 

+ 2^2 q 4 ^ 2cos (g r c) - 2 + qr c sm(qr c )*j (4) 
V(r) = a 2 r 2 : V(p,p',x) = a 2 r 2 5 3 (q) 

+ ^5 (^3gr c cos(gr c ) + (q 2 r 2 c - 3)sin(gr c )^, (5) 

where the potentials are kept fixed at cut-off r c and q — |q| = |p — p'| = \J p 2 + p' 2 — 2pp'x. 
The PW projection of these potentials Vi(p,p') can be obtained by solution of Eq. (2). 
In Tables ^ and ^ we have listed our numerical results for energy eigenvalues. The 
S'-wave energy levels for linear potential are compared with corresponding configuration 
space results and are in excellent agreement with exact energies. As indicated in Table ((T|) 
our numerical results are in excellent agreement with corresponding exact energies E n i = 
(2n + I + ^)hu. In the following we test also the accuracy of our numerical calculations 

TABLE I. Energy eigenvalues for a linear potential. The S*-wave (P and -D-waves) energy levels 
are calculated for m = 1.0 (1.84), ai = 1.0 (0.18263) and r c = 20.0. 



state 


IS 2S 35 45 55 65 


E 
E\3] 

j? _ aiZ t 
(mai ) t 


2.3381 4.0879 5.5205 6.7867 7.9441 9.0226 
2.3373 4.0865 5.5190 6.7814 7.9514 9.0119 
2.3381 4.0879 5.5205 6.7867 7.9441 9.0226 


state 


IP 2P 3P 4P 5P 6P 


E 


0.8830 1.2831 1.6280 1.9454 2.2369 2.5107 


state 


ID 2D 3D AD 5D 6D 


E 


1.1160 1.4789 1.8044 2.1041 2.3844 2.6496 



t Zq are zeros of the Airy function. 



for coulomb potential. PW representation of Fourier transformation of coulomb potential 
V(r) = — a_i/r to momentum space, i.e. V(p,p',x) = j^2, can be obtained analytically 
as Vi(p,p') = Qi( P 2pp> )> wrieT e Qi is the Legendre polynomial of second kind. Clearly 
in the calculation of Vi(p, p') one should overcome the moving singularity which appears in 
Qi at p = p' . To avoid it, one can calculate Vi(p, p') by solution of integral equation (2) 
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and by using the Gauss- Legendre quadrature integration. In Table (Q]) we have compared 
our numerical results for coulomb energy levels with corresponding exact energies. Our 
numerical results confirm the degeneracy of energy levels for different values of I. 

TABLE II. Energy eigenvalues of a quadratic potential for m = 1.0, ai = 0.25 and r c = 10.0. 



state 
E 


15 25 35 45 55 65 75 85 95 
1.500 3.500 5.500 7.500 9.500 11.500 13.500 15.500 17.500 


state 
E 


IP 2P 3P 4P 5P 6P 7P 8P 9P 
2.500 4.500 6.500 8.500 10.500 12.500 14.500 16.500 18.499 


state 
E 


ID 2D 3D 4D 5D 6D 7D 8D 9D 
3.500 5.500 7.500 9.500 11.500 13.500 15.500 17.500 19.499 



TABLE III. Coulomb energy levels for m = 1 and a_i = 1. 



state 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


E 

TP ma -l 

&n — 4rl 2 


-0.2467 
-0.2500 


-0.0619 
-0.0625 


-0.0278 
-0.0278 


-0.1568 
-0.1562 


-0.0101 
-0.0100 



II. HEAVY QUARKONIA MASS SPECTRUM 

In this section we solve the integral equation ([1]) to calculate the mass spectra of heavy 
quarkonia, mesons consisting from heavy quark and antiquark. We consider both linear and 
quadratic confinements. The momentum space representation of the regularized form of 
these potentials can be obtained as: 

V(r) = h air + a : 

r 



V(p, p', x) = I — — + a x r c + a j 5 (q 

-a-i A _ sin(gr c ) \ ^ ai 

2 7r 2 q 2 \ qr c ) 2tt 2 q 

t r / \ ^ — 1 2 

V[r) = h a 2 r + a 



fTTzS 1 - + Trh ( 2 cos (^ rc) - 2 + q r c sin(g r c ) ) , (6) 



r 



V(p,p',x) = (-^r 1 + «2r c 2 + a ^ 3 (q) 



+ Q 9 \ ( 1 - Sm ) + ~~tr~c ( 3g r c cos(g r c ) + (q 2 r 2 c - 3) sin(g r c ) ) . (7) 



2 7r 2 g 2 V qr c J n 2 q 5 



In tables (ITT]) and (ITTJ) the calculated Bottomonium and Charmonium mass spectra are 

n r 

compared with the results obtained by Faustov et al. |3| and also experimental data |4|. 
Our numerical results show that the regularized form of the confining potentials leads to 
energy eigenvalues which are in good agreement with configuration space calculations and 
also experimental data. The study of two-body bound states with other confining potentials 
and also in a relativistic frame is in progress. 

TABLE IV. Charmonium tp(cc) mass spectrum calculated for the sum of linear and quadratic 
confining potentials with the coulomb potential. The parameters of calculation for linear plus 
coulomb (quadratic plus coulomb) potentials are as: ao = —0.29 (—0.05) GeV, a_i = |a s ; a s = 
0.47(0.345), ai = 0.18 GeV 2 (a 2 = 0.174 GeV 3 ), r c = 10.0 (3.0) fm and m c = 1.56 (1.55) GeV . 
State linear+coulomb quadratic+coulomb Exp. [4] 



Bresent Faustov et al. [3] Bresent Faustov et al. [3] 



IS 


3.062 


3.068 


3.076 


3.070 


3.0675 


2S 


3.696 


3.697 


3.720 


3.730 


3.663 


3S 


4.144 


4.144 


4.331 


4.331 


4.159+ 


IP 


3.529 


3.526 


3.492 


3.508 


3.525 


2P 


3.997 


3.993 


4.108 


4.095 




3P 


4.384 


4.383 


4.652 


4.670 




ID 


3.832 


3.829 


3.811 


3.841 


3.770§ 


2D 


4.237 


4.234 


4.396 


4.415 





f 3 Si state 
§ 3 L>i state 
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TABLE V. Same as table (jTIj) but for Bottomonium T(bb) mass spectrum. m& = 4.93 (4.95) GeV 
and a s = 0.39 for linear plus coulomb potentials. Other potential parameters are the same as 
previous table. 

State linear+coulomb quadratic+coulomb Exp. [4] 



Present Faustov et al. [3] Present Faustov et al. [3] 



1 


Q 49^ 


Q 447 


q 7^n 

iy. t OU 


Q 447 


Q 4fif)4t 


2S 


10.006 


10.012 


10.014 


10.007 


10.023^ 


3S 


10.350 


10.353 


10.379 


10.389 


10.355^ 


45 


10.628 


10.629 


10.724 


10.742 


10.580t 


IP 


9.909 


9.900 


9.892 


9.898 


9.900 


2P 


10.263 


10.260 


10.265 


10.259 


10.260 


3P 


10.546 


10.544 


10.594 


10.593 




ID 


10.158 


10.155 


10.135 


10.147 




2D 


10.450 


10.448 


10.488 


10.486 





t 3 5i state 
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